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<3 Abstract 

We discuss the formation of dark matter caustics, and their possible detection by future dark 
^vq matter experiments. The annual modulation expected in the recoil rate measured by a dark matter 

detector is discussed. We consider the example of dark matter particles with a Maxwell-Boltzmann 
velocity distribution modified by a cold stream due to a nearby caustic. It is shown that the effect 
of the caustic flow is potentially detectable, even when the density enhancement due to the caustic 
is small. This makes the annual modulation effect an excellent probe of inner caustics. We also 
show that the phase of the annual modulation at low recoil energies does not constrain the particle 
i - 1 mass unless the velocity distribution of particles in the solar neighborhood is known. 
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I. INTRODUCTION 



Caustics of light have been known since ancient times. The rainbow is a common example 
of a light caustic that forms when the family of refracted light rays is projected on to the 
plane of the sky. Another common example is the heart shaped, or nephroid pattern seen 
on the bottom of a polished coffee cup. Thus caustics of light are regions where the light 
intensity is very large. 
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FIG. 1: Light rays reflected by a polished metal ring or the inner surface of a reflective coffee cup. 
The caustic is the envelope of the family of rays. 

Caustics have well defined physical properties. The light intensity on the inner (concave) 
side of the rainbow varies like the inverse of the square root of distance from the rainbow. 
However, when the rainbow is approached from the outside (convex) side, the light intensity 
remains small until the rainbow is reached, at which point it shoots up abruptly. The coffee 
cup caustic is shown in Fig. 1. The caustic is the envelope of the family of reflected light 
rays, i.e. the curve tangent to all members of the family of rays. There are two qualitatively 
distinct regions separated by the caustic: one with three rays at each point and the other 
with one ray at each point. Consider the point (x — 0,y — 0). Near this point, the light 
intensity varies ~ l^l" 1 when measured along x, and ~ M~ 2 ^ 3 when measured along y. The 
intensity at the point (x = +\e\,y = 0) is double the intensity at the point (x = —\e\,y = 0) 
for small |e|. Near the smooth curve and far from the point (x — 0, y — 0), the light intensity 
only varies in the region with 3 rays at each point, like the inverse square root of the distance 
from the smooth curve. As long as one is far from the boundary between the two regions, a 
small change in position results in only a small change in light intensity. However, close to 
the boundary, a small change in position leads to a very large change in light intensity. 

What we have described are called catastrophes. Caustics are made up of sections of 
catastrophes (more precisely, we are concerned with the bifurcation set of the catastrophe). 
The rainbow caustic is a fold catastrophe. The coffee cup caustic of Fig. 1 consists of 2 fold 
catastrophes that join together forming a cusp catastrophe. Their properties are described 
by catastrophe theory which is a branch of singularity theory applied to physical phenomena. 
Light forms caustics because light is nearly collisionless. Dark matter also has this property. 
We may thus expect dark matter in galactic halos to form caustics. 
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Dark matter caustics are singularities in physical space in the limit of zero velocity dis- 
persion [IH3]. In reality, the finite velocity dispersion of the dark matter particles cuts off 
the divergence in the density. Dark matter caustics are then regions of high density. We 
may expect a continuous infall of dark matter to form caustics provided (i) the particles are 
collisionless and (ii) they have negligible velocity dispersion. 

The shells or arcs seen around some giant elliptical galaxies [U E] are caustics in the 
distribution of stars. When a small galaxy is disrupted, the stars of the small galaxy fall 
into the gravitational potential of the giant galaxy. These stars are then sub-virial, i.e. they 
are cold. Being compact, stars are nearly collisionless. The infall of cold collisionless stars 
leads to the formation of caustics at the turnaround radii. The observations of caustics 
of stars strongly suggests the existence of caustics of dark matter due to the infall of cold, 
collisionless particles. These are called outer caustics and appear as a series of thin spherical 
shells surrounding galaxies. 




(a) (b) 

FIG. 2: Dark matter trajectories, (a) shows purely radial orbits, while (b) shows particles with 
angular momentum. 

Besides the outer caustics, the infall of dark matter particles results in inner caustics. 
When the particles possess angular momentum, they do not pass through the center, and 
hence have non-zero inner turnaround radii. Inner caustics form near these locations. Fig. 2 
shows the dark matter trajectories for the case of (a) radial infall and (b) infall with angular 
momentum. When the trajectories are exactly radial, the infall forms a point singularity at 
the center. When a small amount of angular momentum is added to the flow, the result is a 
caustic made up of sections of catastrophes. Further small perturbations do not change the 
qualitative nature of the inner caustic in (b). The inner caustics are made up of folds and 
cusps which may terminate on higher order catastrophes. The formation of inner caustics 
is described in detail in [6j. Qualitatively, the caustic of Fig. 2(b) resembles the coffee cup 
caustic of Fig. 1. The spherical shell caustics mentioned earlier qualitatively resemble the 
rainbow caustic. 

One may ask whether substructure in real galactic halos can greatly increase the number 
of flows and weaken the caustics. In [7], it was shown that each substructure clump can 
break up a flow into 3 sub-flows, called "daughter flows". Downstream of p clumps, the 
number of sub-flows ~ 3 P can be very large for moderately large p jH]. However this does 
not destroy the discreteness of the phase space hypersurface because the density of the main 
flow is not equally divided among the sub-flows. 

Let us consider the very first flow, i.e. dark matter particles falling towards the halo center 
for the first time. The first turnaround marks the radius at which particles bound to the 
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halo decouple from the Hubble expansion. A caustic does not form at the first turnaround 
radius since there is only one flow there. The first outer caustic occurs at the second outer 
turnaround radius. The first inner caustic occurs at the first inner turnaround radius. We 
thus see that the first inner caustic forms before the first outer caustic. As an example, let us 
assume that the first flow contributes ~ 1% to the local dark matter density in the absence 
of gravitational scattering by substructure. Let us assume there is enough substructure to 
break up the first flow into a large number, say ~ 10 4 of sub-flows. Let us also (incorrectly) 
assume that the density is equally divided among the sub-flows. We then have ~ 10 4 flows, 
each contributing ~ 10~ 4 % to the local dark matter density, and each producing an inner 
caustic of no realistic importance near the inner turnaround radii. The flows eventually fall 
out of the potential well (i.e. their radial velocities becomes positive), and form ~ 10 4 outer 
caustics, of negligible density contrast at the outer turnaround radii. The conclusion then 
would be that there are no physically relevant caustics, inner or outer. 

This conclusion contradicts the observations of P, [10] who have obtained evidence for 
the existence of a dark matter outer caustic in their study of galaxies in the group around 
NGC 5846, and have interpreted the caustic as that of second turnaround, i.e. the first 
outer caustic. This result has been interpreted by [TT] as evidence that an inner caustic 
exists in that group. Indeed it is hard to understand how the first outer caustic can form 
if the infalling flow is so diffused that it cannot produce the first inner caustic. In [3l [8], 
it was found that gravitational scattering of dark matter particles by inhomogeneities does 
not greatly affect the first few flows unless most of the dark matter is composed of very 
massive (~ 1O 1O M ) clumps. We know of no evidence that the dark matter in the Milky 
Way is composed of such massive clumps. On the contrary, the recent findings of p2] show 
that « 40% of the dark matter in halos is accounted for by smooth accretion of particles 
not bound in halos. 

Assuming a continuous infall of cold dark matter particles with angular momentum (as- 
sumed aligned with the baryon angular momentum), one may expect to find a series of inner 
caustics in the plane of the galaxy. Authors [T3HT5] (motivated by the work of |16j [17] ) have 
proposed a model of the Milky Way halo that predicts the locations of the caustics, and the 
approximate density near the caustics. In this model, the earth naturally lies between two 
inner caustics. The distance to the nearest inner caustic has been estimated in [18J. A ring 
like feature in the cluster CI 0024+1654 was studied by [TT] who have interpreted the feature 
as a dark matter inner caustic. The presence of inner caustics in spiral galaxies was inves- 
tigated by [19], who found possible evidence for their existence by analyzing the rotation 
curves of many galaxies. Rises in the rotation curve of the Milky Way also show possible 
evidence for the existence of caustics in our galaxy [18] . In [18] , the presence of a triangular 
feature in the infrared spectrum was proposed as possible evidence for the presence of a 
nearby inner caustic. 

Certain high resolution simulations see cold streams [20J that are associated with inner 
caustics. However, recent N-body simulations performed by [21] do not find significant 
inner caustics. The absence of inner caustics in these simulations is inconsistent with the 
predictions of [HI IT3HT5] and with the claimed observational evidence. In this article, we 
do not attempt to reconcile the simulations with the observations. Instead we suggest that 
a nearby dark matter caustic may be revealed by the velocity distribution of dark matter 
particles in the solar neighborhood. If there exists a nearby dark matter inner caustic, the 
dark matter velocity distribution is influenced by the cold flow forming the caustic. We 
show that the annual modulation expected in the recoil rate measured by a detector can be 
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used to detect the presence of a nearby inner caustic, even for weak density enhancements. 



II. FORMATION OF CAUSTICS 



It was shown in [22J that the continuous infall of cold collisionless matter necessarily pro- 
duces an inner caustic. Let us review the argument: Consider a sphere with a conveniently 
chosen radius, such that all particles of a given flow pass through the sphere. Then, we may 
identify every particle of the flow by 3 parameters (r, 9, <fi). r is the time when the particle 
crossed the sphere on its way in to the halo. 9 and <fi are co-ordinates on the sphere. Let 
x = (x, y, z) be a location in physical space where the dark matter density is measured. A 
catastrophe forms at points where the mapping from the space (r, 9, 4>) to the space (x, y, z) 
is singular. The condition for the mapping to be singular at (x, y, z) is the vanishing of the 
Jacobian determinant: 

Let (9q, (po) be the point on the reference sphere where the angular momentum is a maximum 
(the angular momentum field defined on a sphere has a maximum somewhere unless it is 
zero everywhere). Let us choose xo(to, 6o, <po) to be the point of closest approach of the 

particles with the most angular momentum, and let us define r = \Jxq + y^ + Zq. Then we 
have at Xq. 

dr A dx 

dr dr 

dr A dx 

d9 =X 'd9 = 

dr A dx 

The first condition is due to the reversal of the sign of the radial velocity at xq. The other 
two conditions are due to the angular momentum maximum at (9 ,<f) Q ). Since the four 
vectors dx/dr, dx/d9, dx/d(p, and x cannot all be mutually perpendicular at the same 
point, at least two of the vectors are linearly dependent, or at least one of them is zero. 
This satisfies the caustic condition Eq. [TJ Thus the infall of a cold, collisionless shell of 
particles necessarily results in the formation of an inner caustic, regardless of assumptions 
of symmetry. 

The structure of the inner caustic depends on the spatial distribution of angular mo- 
mentum on the initial reference sphere. For purely curl like initial conditions, the inner 
caustic resembles a ring, made up of sections of the elliptic umbilic catastrophe [6]. For 
purely gradient like initial conditions, the inner caustic has a more complicated shape and 
is made up of sections of the hyperbolic umbilic catastrophe [22]. For the special case of 
axial symmetry with gradient like initial conditions, one may find swallowtail and butterfly 
catastrophes [22] . 

Dark matter caustics can affect the distribution of stars in the galaxy. In [23], it was 
proposed that the Monoceros ring of stars may have formed due to the the presence of the 
n = 2 inner caustic of the Milky Way. Two mechanisms were found by which the star density 
could be enhanced in the vicinity of an inner caustic, and these may have played a role in 
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the formation of the Monoceros ring. Authors [24T - T2T] examined the possibility of observing 
caustics by gravitational lensing, while [2BTI34"] have studied dark matter annihilation in 
caustics. These possibilities however require high densities in the vicinity of the caustic. 

Let us now review the annual modulation expected in the recoil rate of a detector, and 
show that it is an excellent probe of a nearby dark matter caustic even when the density 
enhancement due to caustic is low. It was shown by (35J [36] that due to the motion of the 
earth about the sun, the flux of dark matter particles reaching the earth is modulated with 
a period of 1 year, and is largest at certain times of the year. DAMA is a dark matter direct 
detection experiment that measures the annual modulation in the recoil rate, and claims a 
positive result with 8.9a confidence (37]. Unfortunately, it seems the results obtained by 
DAMA are incompatible with the null results of other experiments (see for example j3BT[4"0~] ) . 

The idea of using the annual modulation as a probe of a nearby dark matter caustic is 
not new. Authors |HH1S] have calculated the expected recoil rate for the late infall self 
similar caustic model. The effect of cold streams has been studied by several authors (see 
for example [46--50J). Here we review the known results, and study the case of a Maxwellian 
velocity distribution modified by a single cold flow due to the presence of a nearby caustic. 
We present results for two cases (i) when the stream contributes significantly to the dark 
matter density (i.e. caustic with a large density boost) and (ii) when the stream is a small 
perturbation. It is shown that even in case (ii), the effects are potentially observable. 



III. DIRECT DETECTION OF DARK MATTER. 

Let us consider Weakly Interacting Massive Particles (WIMPs) as the dark matter. The 
number of recoils per unit time per unit nuclear mass is given by 

dR=—^ [vf(v)dv] , (3) 

tbn m x dQ 

where m-^ = Am p is the mass of the nucleus of the detector, A is the atomic mass number, 
m p is the nucleon mass, p x is the dark matter density at the location of the detector, a is 
the scattering cross section, Q is the recoil energy, and v is the velocity of the dark matter 
particles with respect to the detector. For elastic scattering, the energy transferred to the 
nucleus is given by 

Q = ^(l-cos0), (4) 

where wr = m x m^/(m x + mjsr) is the WIMP-nuleus reduced mass and 9 is the scattering 
angle in the center-of-momentum frame. From Eq. [4], we see that in order to observe a 
nuclear recoil with energy Q, the WIMPs must posses a minimum speed 
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V ^ m K 

since < (1 — cos 9) < 2. It is possible for WIMPs to have a larger speed than the minimum 
given in Eq. [5j but they cannot have a smaller speed for a given recoil energy Q. Here, 
we consider only the spin-independent cross section. The differential cross section da /dQ is 
expressed in the form: 

da (To r,2(r^\ ^O^N r 2 



F 2 (Q) = ^F\Q). (6) 



dQ Qmax 1m\v 2 
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F(Q) is the nuclear form factor. It is obtained from nuclear physics experiments [5TH53] . 
and contains the momentum dependence of the cross section: 



F{Q) = e -U«*)\ (7) 

qr 

in units where h and c are set to 1. q = \/2Qm^^ s = 1 fm, R = 1.2 A 1 / 3 fm, r = \/R 2 — 5s 2 , 
and ji is the spherical bessel function. It is easy to see that F 2 {Q) m 1 for small Q, and 
falls off at large Q. 

We therefore have the number of recoils per unit energy, per unit time, per unit detector 
mass: 

™<!gm H r dv M. ( 8) 

dQ 2m%m x Jv min (Q) v 
Let us define the two quantities: 



POO 

T(Q,t) = / dv 



M(Q 1 ,Q 2 ,t) = / dQ F 2 (Q)T(Q,t). (9) 

T(Q,t) is called the mean inverse speed and M(Qi,Q 2 ,t) is the number of recoils per unit 
time, per unit detector mass, between energies Qi and Q 2 . 

Let us now focus our attention on the distribution of WIMP velocities. The spectrum of 
dark mater velocities is discrete, i.e. a sum over flows: 

/00 (io) 

i 

where Vi is the velocity of the i th flow. There are a large number of flows that occupy the 
inner region of phase space, while the fast flows which are few and well separated in velocity, 
are mostly found in the outer regions of phase space. A detector with finite resolution will 
likely not be able to resolve the numerous flows in inner phase space. It is useful then to break 
up the sum over all flows into 2 parts: a sum over cold flows, and a sum over thermal flows. 
We may replace the sum over thermal flows by a Maxwell-Boltzmann velocity distribution: 

f , 7 _ exp [-{v wh /v ) 2 ] 

7T ' Vq 

where we have ignored the effect of the finite escape velocity. The subscript wh stands for 
"WIMP-halo" and indicates that the velocities are measured in the rest frame of the halo, 



not the detector. We take v to be 220 km/s for the Milky Way. Eq. 10 can be expressed 
as: 

f(v) = S(v-Vi) + / max (t?) = /cold(^) + /max(#)- (12) 

cold flows i 

Let us now compute the mean inverse speed T(Q,t) for the two distributions. Consider 
a single cold flow: f{y w h) = S(v wh — Vfh), where the subscript fh stands for "flow-halo" and 
indicates that the flow velocity Vfh is measured relative to the halo. Transforming to the 
earth's rest frame using Vf h = Vf e + v e (t) + v Q , we obtain f(v we ) = 5{y we — [vf Q — %(£)]). 
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VfQ is the flow velocity relative to the sun, and %>(£) is the velocity of the earth relative to 
the sun. The speed of the flow relative to the earth is given by 



\vfe(t) | = | V{q - % (t) | = v) Q + vl (t) - 2v e (t) ■ v fQ 



1/2 



U /0 



1 - — V @ (t) -Vf Q 
V f@ 



(13) 



for Vq <C f/0. The mean inverse speed is: 

1 



6[\v fe (t)\-v min (Q)\ 



1 



v /0 



fa/0 - % (%(*) ■ *>/©) - ^min(Q)] • (14) 



6{y — vq) is the unit step function = 1 for v > vo and otherwise. t>© is the time averaged 
value of the earth's velocity about the sun. 

Let us now consider the Maxwellian distribution given by Eq. [TTJ We again transform 
to the rest frame of the detector using the relation v w h = v we + v eh = v we + v & + to 
obtain the speed (1-dimensional) distribution f(v e h) = f(v): 



giving a mean inverse speed 



Trv v eh 



( v + v eh \ ' 



(15) 



1 



erf ■ 



2v eh (t) 

with v eh (t) = \v Q + v e (t)\ given by 



^min 

(Q) + v eh (t) 

v 



erf 



Vehif) = vl + V 2 Jt) + 2v e (t) ■ V & 



^min 



1/2 



(16) 



Vq 



H V @ {t) ■ V Q 

Vq 



(17) 



A. Flow near a caustic. 

As we mentioned in the Introduction, the continuous infall of dark matter results in a 
series of outer and inner caustics. If the angular momentum of the dark matter is aligned 
with that of the baryons, these caustics can be expected to lie in the galactic plane. As a 
result, the earth lies between two inner caustics. In the self-similar infall model of [TBHTB] . 
the earth is situated close to the fifth inner caustic. The velocity distribution at the earth's 
location is therefore influenced by the particles forming the caustic. 

Let us choose a co-ordinate system in which the +x axis points towards the galactic 
center, the +y axis points in the direction of galactic rotation, and the +z axis points 
towards the north galactic pole. [15] have listed the velocities of the different dark matter 
flows (note that the system of co-ordinates used here is different from that used in |15]). 
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Here, we will consider only the two flows producing the caustic closest to the sun, which 
have velocities relative to the halo |45] (±100 x + 470 y + z) km/s. Relative to the sun, 
these flows have velocities: 

v aQ = 253.6 km/s (0.3549 x + 0.9345 y - 0.0276 5) 

v bQ = 261.4 km/s (-0.42085 + 0.9067 y - 0.0268 z) . (18) 




Jan 21 May 5 July 22 Nov 4 



FIG. 3: Speed of the two caustic flows relative to the earth. Flow (a) peaks in November, and is 
smallest in May. Flow (b) peaks in January, and is smallest in July. 

The velocities of the sun (about the halo center) and the earth (about the sun) in these 
co-ordinates are respectively (see for example [Ml H7J EO] , and references therein): 

v Q = 233.3 km/s [0.0429 £ + 0.9986 y + 0.0300 ,2] 
= 29.8 km/s [(0.9931 cos - 0.0670 sin 0) x 

+ ( 0.1170 cos + 0.4927 sin 0) y - (0.0103 cos + 0.8676 sin <f>) z] , (19) 

where the angle 0(t) = 2n (t — March 21)/365. We note that earth's velocity about the sun 
is most closely aligned with the sun's velocity about the halo center when = 71°, which 
occurs around June 1. The two velocity vectors are most misaligned when = 251° which 
occurs six months later, around Nov 30. F ig. |3| sh ows the two caustic flow velocities v a ®(t), 



and fbe(t) relative to the earth. Using Eq. 13, 18, and 19, we see that v a ®(t) is largest when 



= 225° (around Nov 4) and smallest when = 45° (around May 5). is largest when 

= 302° (around Jan 21) and smallest when = 122° (around July 22). 

One of these flows is the dominant flow, and for definiteness, let us assume flow v a is the 
dominant flow. This flow shows an w 8.3% speed modulation with a maximum of ~ 275 
km/s on ~ November 4, and a minimum of ~ 232 km/s on ~ May 5. 

Let £ be the fraction contributed by the dominant caustic flow to the total dark matter 
density at the earth's location. Then, 

Pf(v) = P [UnoM + (1 - • (20) 
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£ can be large if the density enhancement due to the caustic is large. [HJ HH] have estimated 
that the dominant caustic flow, dubbed the "big flow" contributes ~ 73% of the dark matter 
density at the earth's location. We now show that even if the dominant flow contributes as 
little as 5% to the local dark matter density, it has potentially observable consequences. 
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FIG. 4: Mean inverse speed T(Q,t). Panel (a) shows the variation for the caustic flow (a) and for 
the Maxwellian. Panel (b) shows a combination of the two. Shown are measurements in November 
and June. 



B. Annual modulation. 



Fig. [4j shows the mean inverse speed T(Q,t) for an assumed WIMP mass m x = 70 GeV, 
and a Germanium detector with atomic mass number A = 73, as a function of recoil energy 
Q. The top panel (a) shows T(Q,t) for the two extreme cases £ = 1 (density entirely due 
to the caustic flow, shown by the two solid curves black and red) and £ = (caustic flow 
absent, shown by the dashed curves blue and pink), expected in November, and in June 
respectively. For a flow velocity relative to the earth vf, the maximum energy at which 
recoils are measured is given by 

2/7? ^ 1)^ 

Qmax = — — « 25keV [1 ± 0.166] , (21) 
m N 

for the assumed values of m x , m^, and Vp, with the maximum and minimum values occurring 



in November and May respectively. The height of the step is given by Eq. 14 and peaks 
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when the speed is the lowest i.e. in May, and is smallest in November when the speed is 
greatest. Note that the modulation of the edge of the step (oc square of the flow speed) is 
twice the modulation of the step height (inverse of the flow speed), and has opposite phase. 
The flow ceases to be visible to the detector during part of the year starting from Q ~ 21 
keV (for the assumed values of m x and m^). As the energy is increased, T is non-zero close 
to November, resulting in a maximum then. 

Let us look at the behavior of the recoil rate integrated over all energies M(0, oo, t): 



M(0,oo,t) = / dQF 2 (Q)T(Q,t) = dQ 
Jo Jo 



F 2 (Q) 



Qmax /r>r>\ 
~ VfQ, (22) 

V f@ 

where we ignored the energy dependence of F 2 (Q). Thus, when integrated over all recoil 
energies, the recoil rate peaks in November for this particular flow, and is smallest in June. 
However, when M is measured over only small recoil energies, we have: 

M(0, Q small , t) = / Qsma " dQF\Q) T(Q, t) ~ ^ mn ^, (23) 
Jo v fQ 

which implies that the integrated recoil rate peaks in June for small recoil energies Q. As 
Q is increased, the phase abruptly reverses, peaking in November. The abrupt shift in T is 
because the flow is no longer visible to the detector around June at these energies, owing 
to the smaller speed around June. In contrast, for a Maxwellian velocity distribution, the 
recoil rate peaks in November for small Q, reverses phase as Q is increased, and peaks in 
June for sufficiently large Q [53] . 

The two dashed curves show the modulation of T for the purely Maxwellian case. The 
difference is far less dramatic than for the case of the pure caustic flow. At low recoil energies, 
T peaks in November (blue dashed curve) and is smallest in June (pink dot-dashed curve). 
The phase reverses at ~ 15 keV (for the assumed values of m x = 70 GeV and A = 73), and 
T then peaks in June and has a minimum in November. At high energies, T is small since 
(i) there are few particles with sufficient energy to effect recoils and (ii) F 2 (Q) is small for 
large Q. 

Fig. |4](b) shows the more physical case with both the caustic flow and the thermal halo 
contributing to the recoil energy spectrum. Shown are the two cases £ = 0.5 (strong density 
enhancement) and £ = 0.05 (weak density enhancement). Let us first consider the case 
£ = 0.5 (shown by the solid red and black curves). The caustic flow dominates until the 
flow ceases to be visible to the detector during parts of the year, resulting in a sharp drop in 
T. For the case of £ = 0.05, the difference in T between June and November measurements 
is not so large. At low energies, T is mostly due to the Maxwellian, and is largest in 
November. With increase in recoil energy, the modulation amplitude decreases, reversing 
phase at ~ 12.5 keV. The peak now occurs in June. The next change occurs at Q ~ 22 
keV, with T measured in November being larger than in June, owing to the absence of the 
caustic flow near June. The phase changes again at Q ~ 29 keV, due to the total absence 
of the caustic flow. For Q > 29 keV, the modulation is consistent with a Maxwellian halo, 
with a maximum measured in June. 

Fig. [5] shows the annual modulation in M(Q±, Q2, t) as a function of the time of year, for 
the case £ = 0.5 (caustic flow contributes 50% to the local dark matter density). The six 
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panels show different recoil energy ranges, with a bin size of 5 keV. Panels (a), (b), and (c) 
are qualitatively identical, and the caustic flow is dominant, with the recoil rate being largest 
in June, and smallest in November. The caustic flow is not seen during part of the year in 
Panel (d), resulting in the steep minimum around June. The modulation is about 25%. In 
panel (e), much of the flow is absent around June, leading to a near ~ 50% modulation in 
recoil energy. The flow is completely invisible in panel (f), and the modulation of ~ 3% is 
due to the Maxwellian distribution. 




Jan Apr Jul Oct Dec Jan Apr Jul Oct 



FIG. 5: Annual modulation in M(Qi,Q2,t) for the caustic flow (a). Q is in keV. Case £ = 0.5. 

We see from panel (a) that M peaks in June and is lowest in November, in the 5 — 10 keV 
energy band. If one were to incorrectly assume a Maxwellian distribution with no cold flow 
contribution, panel (a) would imply an upper limit for the WIMP mass of about m x = 30 
GeV (for the assumed A = 73), far below the value of m x = 70 GeV used in this example. 
Thus, the phase of the annual modulation constrains the WIMP mass [54J only when the 
distribution of particle velocities in the solar neighborhood is known. One may however use 
the results of two different experiments to constrain the mass without assuming a form for 
the velocity distribution jSHJ EH] . 

Fig. [6]shows M(Qi, Q 2 , t) for the case £ = 0.05 (a caustic flow contribution of 5%). In the 
5—10 keV energy range, M shows a minimum in June, and a maximum in November due 
to the dominant contribution of the Maxwellian halo. Phase reversal occurs in the 10 — 15 
keV region, and the modulation is smallest at these energies. In panel (c), M peaks in June 
and is minimum in November. The caustic flow influences M in panels (d) and (e) because 
at these energies, the flow is visible to the detector only during certain months of the year. 
As a result, even a ~ 5% contribution from the caustic flow can lead to a ~ 5% modulation 
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effect, comparable to the modulation produced by the dominant Maxwellian component. In 
panel (e), the peak of M occurs in November, and the modulation is far from sinusoidal. 
In panel (f), the caustic flow does not contribute, and the recoil rate is purely due to the 
Maxwellian distribution. 




Jan Apr Jul Oct Dec Jan Apr Jul Oct 



FIG. 6: Annual modulation in M(Qi,Q2,t) for the caustic flow (a). Q is in keV. Case £ = 0.05. 



IV. CONCLUSIONS 

We have discussed the formation of dark matter caustics and provided possible observa- 
tional evidence for the existence of outer and inner caustics. The self-similar infall model 
predicts the locations of the inner caustics, and their approximate densities. According to 
this model and the observations, the earth may be located near the fifth inner caustic. One 
may hope to detect the nearby caustic by means of the annual modulation expected in the 
recoil rate of a dark matter detector. 

We discussed the annual modulation effect and derived expressions for the expected recoil 
rate. We considered the case of a Maxwellian velocity distribution modified by a cold flow 
due to the nearby caustic in the self-similar infall model. Fig. 3 shows the velocities of 
the two flows that contribute to the formation of the nearest caustic. Fig. 4 shows the 
mean inverse speed T(Q,t) expected for a cold flow, for a Maxwellian distribution, and for 
a combination of the two. We then computed the recoil rate integrated over energies in 5 
keV bins. We considered two cases: (i) a 50% contribution due to the caustic flow (Fig. 5) 
and (ii) a 5% contribution (Fig 6). We showed that the phase of the annual modulation 
at low energies can be used to constrain the WIMP mass only if the velocity distribution 
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of dark matter particles in the solar neighborhood is known. We also showed that even a 
small contribution by a caustic flow can significantly alter the modulation of the recoil rate, 
at energies near the threshold energy of the caustic flow. Annual modulation is thus an 
excellent tool to detect a nearby dark matter caustic. 



Acknowledgments 

A.N. acknowledges financial support from the Bruce and Astrid Mc Williams postdoctoral 
fellowship. 



[1] V.I. Arnold, S.F. Shandarin, Y.B. Zeldovich, Geophysical and Astrophysical Fluid Dynamics, 
vol.20, p.lll (1982) 

S.F. Shandarin, Y.B. Zeldovich, Reviews of Modern Physics, vol.61, p. 185 (1989) 
P. Sikivie, J.R. Ipser, Physics Letters B, vol.291, p. 288 (1992) 

D. F. Malin, D. Carter, Nature (London), vol.285, p. 643 (1980) 
L. Hernquist, P.J. Quinn, The Astrophysical Journal, vol.312, p.l (1987) 
P. Sikivie, Physical Review D, vol.60, p.063501 (1999) 
P. Sikivie, S. Wick, Physical Review D, vol.66, p.023504 (2002) 
A. Natarajan, P. Sikivie, Physical Review D, vol.72, p. 083513 (2005) 
A. Mahdavi, N. Trentham, B.R. Tully, The Astronomical Journal, vol.130, p. 1502 (2005) 
N. Trentham, B.R. Tully, A. Mahdavi, Monthly Notices of the Royal Astronomical Society, 
vol.369, p.1375 (2006) 

V. Onemli, P. Sikivie, Physics Letters B, vol.675, p. 279 (2009) 

S. Genel, N. Bouche, T. Naab, A. Sternberg, R. Genzel, |arXiv:1005.4 058v2 [astro-ph] 
P. Sikivie, I.I. Tkachev, Y. Wang, Physical Review Letters, vol.75, p.2911 (1995) 
P. Sikivie, I.I. Tkachev, Y. Wang, Physical Review D, vol.56, p.1863 (1997) 
P. Sikivie, L.D. Duffy, Physical Review D, vol.78, p. 063508 (2008) 
J. A. Filmore, P. Goldreich, The Astrophysical Journal, vol.281, p.l (1984) 

E. Bertschinger, The Astrophysical Journal Supplement Series, vol.58, p. 39 (1985) 
P. Sikivie, Physics Letters B, vol.567, p.l (2003) 
W.H. Kinney, P. Sikivie, Physical Review D, vol.61, p. 087305 (2000) 
D. Stiff, L.M. Widrow, Physical Review Letters, vol.90, p.211301 (2003). 
M. Vogelsberger, S.D.M. White, R. Mohayaee, V. Springel, Monthly Notices of the Royal 
Astronomical Society, vol.400, p.2174 (2009) 

A. Natarajan, P. Sikivie, Physical Review D, vol.73, p. 023510 (2006) 
A. Natarajan, P. Sikivie, Physical Review D, vol.76, p. 023505 (2007) 
C. Hogan, The Astrophysical Journal, vol.527, p. 42 (1999) 

C. Charmousis, V. Onemli, Z. Qiu, P. Sikivie, Physical Review D, vol.67, p. 103502 (2003) 
R. Gavazzi, R. Mohayaee, B. Fort, Astronomy and Astrophysics, vol.445, p. 43 (2006) 
V. K. Onemli, Physical Review D, vol.74, 123010 (2006) 
C. Hogan, Physical Review D, vol.64, p.063515, (2001) 

R. Mohayaee, S. Shandarin, Monthly Notices of the Royal Astronomical Society, vol.366, 
p.1217 (2006) 



14 



[30] R. Mohayaee, S. Shandarin, J. Silk, Journal of Cosmology and Astroparticle Physics, Issue 
05, p.15 (2007) 

A. Natarajan, Physical Review D, vol.75, p. 123514 (2007) 
A. Natarajan, P. Sikivie, Physical Review D, vol.77, p. 043531 (2008) 
N. Afshordi, R. Mohayaee, E. Bertschinger, Physical Review D, vol.79, p. 083526 (2009) 
S.D.M. White, M. Vogelsberger, Monthly Notices of the Royal Astronomical Society, vol.392, 
p.281 (2009) 

A.K. Drukier, K. Freese, D.N. Spergel, Physical Review D, vol.33, p. 3495 (1986) 
K. Freese, J. Frieman, A. Gould, Physical Review D, vol.37, p.3388 (1988) 
R. Bernabei et al. for the DAMA collaboration, The European Physical Journal C, vol.67, p. 
39 (2010) 

M. Fairbairn, T. Schwetz, Journal of Cosmology and Astroparticle Physics, Issue 01, p. 037 
(2009) 

C. Savage, K. Freese, P. Gondolo, D. Spolyar, Journal of Cosmology and Astroparticle Physics, 
Issue 09, p. 036 (2009) 

C. Savage, G. Gelmini, P. Gondolo, K. Freese, |arXiv:1006.0972"l [astro-ph] (2010) 

C. J. Copi, L.M. Krauss, Physical Review D, vol.63, p.043507 (2001) 
A.M. Green, Physical Review D, vol.63, p. 103003 (2001) 
J.D. Vergados, Physical Review D, vol.63, p. 063511 (2001) 
G. Gelmini, P. Gondolo, Physical Review D, vol.64, p.023504 (2001) 

F. S. Ling, P. Sikivie, S. Wick, Physical Review D, vol.70, p. 123503 (2004) 

D. Stiff, L. Widrow, J. Frieman, Physical Review D, vol.64, p.083516 (2001) 
K. Freese, P. Gondolo, H. Newberg, M. Lewis, Physical Review Letters, vol.92, p. 11 1301 
(2004) 

P. Gondolo, G. Gelmini, Physical Review D, vol.71, p. 123520 (2005) 
K. Freese, P. Gondolo, H.J. Newberg, Physical Review D, vol.71, p. 043516 (2005) 
C. Savage, K. Freese, P. Gondolo, Physical Review D, vol.74, p.043531 (2006) 
R.H. Helm, Physical Review, vol.104, p. 1466 (1956) 
J.D. Lewin, P.F. Smith, Astroparticle Physics, vol.6, p. 87 (1996) 

G. Jungman, M. Kamionkowski, K. Griest, Physics Reports vol.267, p. 195 (1996) 
M.J. Lewis, K. Freese, Physical Review D, vol.70, p.043501 (2004) 

M. Drees, C.L. Shan, Journal of Cosmology and Astroparticle Physics, Issue 06, p. 011 (2007) 
M. Drees, C.L. Shan, Journal of Cosmology and Astroparticle Physics, Issue 06, p. 012 (2008) 



15 



